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Abstract 


Mean  flow  and  compressible  turbulence  measurements  have  been  obtained  at  two 
stations  downstream  of  low-angle  supersonic  injection  into  a  Mach  3.0  flow  (Re/m  =  15  x 
10*’).  Data  were  collected  using  conventional  Pitot  and  cone-static  probes,  single  and 
multiple  overheat  cross-wire  anemometry,  and  flow  visualization  techniques  (shadowgraphs 
and  schlierens).  A  direct  measure  of  total  Reynolds  shear  stress  was  accomplished  using  a 
turbulence  transformation  of  the  Reynolds-averaged  Navier-Stokes  equations.  Results  show 
compressibility  effects,  indicated  by  density  fluctuations,  to  be  large  relative  to  the  velocity 
fluctuations  and  on  the  same  order  in  all  three  components.  Compressibility  appears  to 
account  for  as  much  as  75%  of  the  total  Reynolds  shear.  The  results  of  the  present  study 
suggest  turbulent  compressibility  effects  are  very  important  for  this  class  of  flows. 


Compressible  Turbulence  Measurement  in  Low-Angle 
Injection  into  a  Supersonic  Flow 


L  Background 


1.1  Introduction 

One  of  the  linchpins  of  the  once  (and  future?)  National  AeroSpace  Plane  is  its 
supersonic  combustion  ramjet  (SCRAMJet)  engines.  These  engines  require  fuel  to  be 
injected  into  a  supersonic  flow,  rapidly  mixed,  and  burned.  Because  of  the  expense  and 
difficulty  in  operating  such  an  engine,  accurate  computer  models  of  the  engine  and 
associated  flow  physics  are  required  to  optimize  the  design. 

Some  types  of  high-energy  gas  lasers  contain  flowflelds  similar  to  that  of  a 
SCRAMjet  engine.  One  such  laser  is  the  chemical  oxygen-iodine  laser  (COIL)  being 
developed  as  the  heart  of  the  Airborne  Laser  (ABL)  theater  missile  defense  platform  under 
development  by  the  Air  Force. 

Whether  a  SCRAMjet  engine  or  a  high-power  gas  laser,  the  complex  flowfield  within 
these  devices  must  be  understood  for  the  concepts  to  reach  fruition.  Building  models  and 
mockups  and  testing  them  in  wind  tunnels  is  time  consuming,  expensive,  and,  for  the 
extreme  temperature  cases,  not  feasible.  Thus,  researchers  and  engineers  are  relying  more 
heavily  on  computational  techniques  to  develop  and  optimize  these  concepts.  However, 
these  techniques  are  hindered  by  inaccurate  turbulence  models  which  require  accurate 
compressible  turbulence  data  for  verification  and  improvement.  Presently,  accurate 
compressible  turbulence  data  is  scarce. 

1.2  Requirements  for  Turbulence  Modeling 

Large  supersonic,  high  temperature  wind  timnels  are  expensive  and  have,  in  the  last 
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thirty  years,  become  more  scarce.  Since  1963,  fiilly  two-thirds  of  these  types  of  wind  tunnels 
have  been  dismantled  or  abandoned.'  At  the  same  time,  the  cost  of  computing  power  has 
decreased  several  orders  of  magnitude.  But  computers  are  not  a  panacea.  Computers,  it 
should  be  noted,  do  exactly  what  they  are  programmed  to  do:  a  code  that  executes 
successfully  does  not  necessarily  produce  accurate  answers. 

In  the  field  of  high  speed  turbulence  modeling,  the  complete  solution  of  the  unsteady 
Navier-Stokes  equations  is  the  ultimate  goal.  Today’s  computers,  though  powerful,  still  lack 
the  speed  and  memory  required  for  such  a  complete  solution.  Although  estimates  vary,  it 
is  commonly  held  that  10  grid  points  are  necessary  for  adequate  resolution  of  the  motion  of 
a  turbulent  eddy.^  The  relative  size  of  large  scale  and  small  scale  flow  structures  is 
commonly  given  by' 

Large  scale  _  ^^3/4 
Small  scale  ^ 

For  the  current  work.  Re/m  =  15x  lO'^^  and  6  ~  1  cm  .-.  Re6=  15x10^.  The  ratio  between  size 
of  large  and  small  scale  structures  for  the  present  work  is  7600:1.  This  would  require 
approximately  76,000  points  across  the  boundary  layer  and  (76,000)'  points  for  1  cnf  of 
flow.  This  enormous  grid  requirement  must  be  coupled  with  an  equally  challenging 
temporal  resolution.  It  is  unlikely  the  direct  numerical  simulation  of  a  high  Reynolds 
number,  compressible,  turbulent  flow  field  will  yield  completely  to  computational  methods 
in  the  near  future. 

Engineering  efforts  in  the  near  term  will  have  to  focus  on  accurate  modeling  of  the 
small-scale  compressible  turbulence.  The  two  primary  avenues  to  this  end  are  the  Reynolds 
(time)-averaged  Navier-Stokes  equations  (RANS)  and  the  Favre  (mass-weighted)-averaged 
Navier-Stokes  equations  (FANS)  (both  will  be  discussed  further  in  Chapter  2).  Both 
methods  require  closure  formulations. 

The  development  of  the  needed  closure  formulations  (turbulence  models)  rely 
heavily  on  both  experimental  data  and  physical  intuition.  The  lack  of  compressible 
turbulence  data’*  has  lead  to  direct  extensions  of  incompressible  models  and  sometimes  ad- 
hoc  assumptions.  These  extensions  and  assumptions  have  generally  shown  to  be  invalid  for 
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compressible  flows. ^  The  availability  of  experimental  data  has  been  hampered  by  the 
expense  and  difficulty  of  measuring  the  compressible  components  of  the  turbulence.  In 
addition,  it  is  necessary  to  employ  experimental  methods  which  can  measure  the  needed 
compressible  terms  in  the  Navier-Stokes  equations  and  to  report  the  data  in  a  form 
applicable  to  both  theoretical  and  computational  efforts.  These  issues  have  been  addressed 
by  Bowersox®  and  are  employed  in  the  present  work. 

1.3  Scope  of  Present  Study 

The  present  study  seeks  to  expand  upon  the  current  understanding  of  mixing  and 
turbulence  in  a  low-angle  supersonic  injection  into  a  supersonic  flow.  A  simple  injection 
model,  based  on  parameters  of  similar  experiments  in  the  literature  is  employed.  Mean 
flow  and  single  overheat  cross-wire  anemometry  is  taken  at  two  downstream  locations  to 
provide  data  which  will  lead  to  a  more  accurate  determination  of  relative  importance  of  the 
terms  in  the  RANS.  Multiple  overheat  cross-wire  anemometry  is  taken  at  a  single  location 
to  validate  assumptions  made  in  the  taking  of  single  overheat  data.  The  ultimate  aim  of  the 
present  work  is  to  provide  data  which  will  aid  the  validation  of  CFD  codes,  improve 
compressible  turbulence  modeling,  and  enhance  the  understanding  of  the  associated  flow 
physics. 


IL  Mathematical  Development 

This  chapter  presents  the  development  of  the  forms  of  the  unsteady  Navier-Stokes 
equations  utilized  in  the  present  work.  As  stated  in  Chapter  1,  the  two  most  common  forms 
implemented  are  the  Reynolds  (time)-averaged  (RANS)  and  the  Favre  (mass-weighted)- 
averaged  Navier-Stokes  (FANS)  equations. 

2. 1  Compressible  Navier-Stokes  Equations 

The  full  unsteady,  compressible  Navier-Stokes  equations^  are  presented  below  as  a 
prelude  the  development  of  the  RANS  and  FANS; 


■^-£-  +  -^(p  u]  =  0 
dt  dx/ 

(2.1) 

a,  ^  a,  .  dP  ^  d  .. 

(2.2) 

(2.3) 

where 


K  =  +  I ",  "/ 


T..  =  U 

y  ^ 


d  u.  d  Uj 

V  d-v 


+  6..X 


d  u. 


''  ax. 


q.  =  -k 


dT 

dx, 


(2.4) 

(2.5) 

(2.6) 
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(2.7) 


A  =  -  —  |a  {hulk  viscosity) 
3 


As  stated  in  Chapter  1,  numerical  solutions  to  the  above  equations  are  not  available 
at  this  time,  therefore  approximate  methods  must  be  utilized.  Time  averaging  of  the  above 
equations  leads  to  the  RANS  equations. 


2.2  Reynolds- Aver  aged  Navier-Stokes  (RANS)  Equations 

A  standard  method  of  dealing  with  turbulent  flows  is  to  separate  the  instantaneous 
flow  properties  into  mean  and  fluctuating  components.  The  instantaneous  flow  properties 
become: 


u.  =  u.  +  u. 

I  I  I 

P^F^P' 

p  =p  +  p' 


<7, 


(2.8) 


T.. 

'/ 


T..  +  T,. 

'/  '/ 


The  barred  quantities  are  time  averaged  over  a  period  longer  than  the  longest  fluctuation  but 
shorter  than  any  unsteadiness  of  the  mean  flow.  The  prime  quantities  are  the  deviations 
from  the  average. The  flow  variables  are  time  averaged  according  to 


'0 

w=  I  u{x,y,z. 


t)dt 


(2.9) 


By  substituting  the  mean  and  fluctuating  components  (2.7)  into  the  governing 
equations  (2.  l)-(2.3)  and  employing  (2.9),  the  RANS  equations  are  developed 


iP 

dt 


+  -^(p  u  +  p' w7)  =  0 

j 


(2.10) 
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(2.11) 


A(p  „  +  p  '  uf)  +  ^(p  u.  u)  = 
dt  ^  ^  dXj 


dP 

d  X. 


d  t 


(P 


9'K') 


^  (p/^o",)  = 


d  X 


dx 


+  U.  T.. 
y  '  y 


(2.12) 


The  compressible  terms  in  (2.  l0)-(2. 12)  can  be  written  as 

m.  = -o'  u/ 

I  r  I 


(2.13) 


=  -  P  w  '  w  '  -  U  p'  -  U- p'  M.'  -  p'  u  '  u' 

y  I  ^  J  J  ‘  ^  I  J 


(2.14) 


T 


+  /Zq  P'  w,'  +  W  P'  /^o'  +  p'  K' 


(2.15) 


where  m  ^  is  the  turbulent  apparent  mass,  is  the  turbulent  compressible  shear,  and  q  ^ 
is  the  turbulent  compressible  heat  flux.  The  averaging  process  is  notable  in  that  it  removes 
some  flow  information,  in  particular,  its  frequency,  phase,  and  wavelength. 

In  incompressible  flow,  the  apparent  mass  and  the  last  three  terms  of  both  (2. 14)  and 
(2.15)  are  identically  zero.  In  compressible  flows,  the  last  term  of  both  (2. 14)  and  (2.15) 
is  a  triple  correlation  and  therefore  assumed  small  and  negligible.  The  second  and  third 
terms  are  significant,  however,  and  their  exclusion  in  incompressible  extensions  has  led  to 
the  latter's  failure  to  fully  model  the  physics  of  compressible  turbulence. 

2.3  Favre- Aver  aged  Navier-Stokes  Equations 

The  second  averaging  method  is  the  Favre-averaged  (mass-weighted)  Navier-Stokes 
equations.  A  thorough  treatise  on  the  development  of  the  FANS  can  be  found  in  Schetz. " 
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Like  the  RANS,  FANS  eliminates  some  flow  information  in  addition  to  the  phase,  frequency, 
and  wavelength,  in  this  case,  the  explicit  moments  of  the  density  fluctuation.  The  FANS  has 
been  applied  to  compressible  flows,  using  extensions  of  incompressible  turbulence  models, 
with  some  success  with  wall  shear  layers,  but  generally  with  failure  with  free  shear  layers. 
The  flow  conditions  in  the  present  work  include  both  wall  shear  and  free  shear,  but  with 
emphasis  on  the  free  shear. 

2.4  Turbulence  Transformation 

As  previously  stated,  different  formulations  of  the  Navier-Stokes  equations  lend 
themselves  to  different  purposes.  Thus,  right-hand  sides  of  equations  (2. 1 3)-(2. 1 5)  must  be 
transformed  into  experimentally  measurable  terms.  The  turbulence  transformation  below 
follows  from  that  developed  by  Bowersox.'^ 

The  following  identity  begins  the  transformation 


p4>  "(p4>)  +  (p4>)' 

(2.16) 

p4>  =(p  +  p')($  +  4>') 

(2.17) 

where  (|)  is  a  generic  variable.  Setting  the  right-hand  sides  of  (2. 1 6)  and  (2. 1 7)  equal  yields 

pcf)  +  (pcj))' =(p  +  p')((p +4)')  (2.18) 

Time  averaging  (2. 17)  yields 

pcj)  =  p(p  +  p'(()'  (2.19) 

which  is  also  a  preliminary  step  in  the  FANS  development.  Subtracting  (2. 1 9)  from  (2.18) 
produces  the  fluctuation  relation 

(p  4))'  =  p  (f)'  +  ^  p'  +  (p'  ({)'  -  p'  (|)') .  (2  20) 

Substituting  4>  =  w,  and  4)  =  Uj  into  (2.20)  and  forming  the  products  (p  w,)’  and  (p  uf  one  can 
collect  the  terms  in  the  Reynolds  shear  (2. 14)  on  one  side  of  the  equation  to  produce  the 
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Reynolds  shear  transformation 


(P  «, )'  (P  «, )' 


QU.U.  — 
^  ‘  j  -z: 


(2.21) 


where  the  triple  correlation  p '  u. '  u. '  and  fourth  order  terms  have  been  neglected.  Similarly, 
the  Reynolds  heat  flux  terms  can  be  shown  to  be 

(2.22) 


All  the  terms  appearing  in  (2.21)  and  (2.22)  can  be  measured  directly  or  estimated  by 
separation  of  cross-wire  variables.  Chapter  4  contains  a  more  detailed  description  of  the 
process. 
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IIL  Facilities  and  Instrumentation 


This  chapter  identifies  the  facilities  and  instrumentation  utilized  in  the  present  study. 
The  AFIT  Maeh  3  wind  tunnel  and  injector  model  will  be  considered  as  will  the  flow 
visualization  techniques  and  equipment.  Data  acquisition  and  instrumentation  for  both  mean 
flow  and  turbulence  (hot-wire)  will  also  be  covered. 

3. 1  Supersonic  Wind  Tunnel 

The  present  experimental  work  was  accomplished  in  the  AFIT  Mach  3  wind  tunnel. 
The  wind  tunnel  has  a  6.35  cm  x  6.35  cm  test  section  and  operates  with  a  combination  of 
pressure  and  vacuum.  The  tunnel  has  a  replaceable  Mach  3.0  nozzle  section  and  a  variable 
diffuser.  The  pressure  side  of  the  tunnel  is  a  0.69  MPa  house  air  supply  provided  by  two 
Atlas  Copco  GAU  807  compressors  delivering  0.5  kg/sec.  Two  Pioneer  Refrigerant  Air 
Dryers  dry  the  house  air  supply.  The  air  is  further  dried  and  cleaned  prior  to  its  introduction 
into  the  tunnel  via  a  cyclone  separator  and  multiple  layers  of  a  eloth  fiber  filter  paper. 

The  vacuum  side  of  the  wind  tunnel  is  provided  by  sixteen  storage  tanks  (~  16  m^ 
total  volume)  evacuated  by  three  7.5  hp,  230  V  Stokes  MicroVac  pumps.  This  arrangement 
could  provide  a  vacuum  of  approximately  2  mmHg,  but  runs  were  typically  made  at  around 
10  mmHg.  The  wind  tunnel  had  a  useful  run  time  of  approximately  20-25  seconds. 

The  tunnel  settling  chamber  contains  a  multi-screen  flow  straightener  with  a 
thermocouple  and  Pitot  probe  downstream  just  prior  to  the  convergent  portion  of  the  nozzle 
to  account  for  any  pressure  drop  across  the  straightener.  An  Endevco  0-690  kPag  pressure 
transducer  monitored  chamber  stagnation  pressure  through  the  Pitot  probe.  The 
thermoeouple  was  used  to  measure  total  temperature  of  the  flow.  An  Omega  Engineering 
Type  K  (chromel-alumel)  constructed  of  Type  316  stainless  steel  was  employed. 

3.2  Injector  Model 

As  stated  in  Chapter  1,  the  purpose  of  the  present  study  was  to  obtain  three- 
dimensional  turbulence  for  a  supersonic,  low  angle  flow  injected  into  a  supersonic 
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FIGURE  3.1;  Tunnel  Coordinates 

freestream  (see  Figure  3.1). 

3.2.1  Injector  Model  Considerations 

There  are  several  critical  parameters  to  be  considered  for  the  injection  model: 
injection  angle,  yaw  angle,  effective  back  pressure,  and  injector  Mach  number.  For  the 
present  study,  the  injector  parameters  were  also  scaled  to  be  compatible  with  existing 
equipment: 


injection  angle  (0) 

25° 

yaw  angle  (p) 

0° 

effective  back  pressure  ratio  (P/Peb) 

3.25 

injector  Mach  number  (Mj) 

1.8 

3.2.2  Injector  Model  Design 

The  injector  (Figure  3.2)  was 
designed  with  ease  of  construction  in  mind. 
The  upper  portion  of  the  air  inlet  is  threaded 
to  accept  a  6.35  mm  (inner  diameter)  pipe 
and  expands  conically  to  19.05  mm 
diameter.  The  settling  chamber  is  a  19.05 
mm  diameter  spherical  bottom  cavity.  A 
hole  was  tapped  into  the  side  of  the  settling 


chamber  to  accept  an  Endevco  0-690  kPag 


FIGURE  3.2;  Injector  Model 
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chamber  to  accept  an  Endevco  0-690  kPag  pressure  transducer.  The  lower  portion  of  the 
injector  block  was  machined  to  be  perfectly  flush  with  the  ceiling  of  the  test  section  to 
reduce  edge  shocks.  The  injector  is  conical  with  a  throat  diameter  of  3.264  mm  and  an 
effective  exit  diameter  (d)  of  3.86 1  mm.  The  dimensions  of  the  injector  exit  are  4.0  mm 
X  9.5  mm. 


3.2.3  Effective  Back  Pressure 

The  concept  of  “effective  back  pressure'XPcb)  was  first  put  forth  by  Schetz  and 
Billig.''^  The  initial  model  was  based  on  treating  the  normally  injected  supersonic  stream 
as  a  cylindrical  object  in  the  freestream.  The  PebWas  one-half  the  stagnation  pressure  on  the 
surface  of  the  cylinder  based  on  the  freestream  conditions.  A  later  work  by  Billig,  et  al. 
set  P^b  at  two-thirds  the  freestream  Pitot  pressure.  Additional  work  set  P„b  equal  to  80%  of 
the  static  pressure  behind  a  normal  shock  in  the  freestream.  There  is  little  difference  in 
the  Peb  obtained  for  these  different  methods.  More  recent  works’^  '**  extended  the  models 
to  angled  injection  (as  opposed  to  the  transverse  injection  of  earlier  works). 

The  new  model  assumes  the  low  angle  jet  to  behave  like  a  cone  whose  half  angle  is 
the  angle  of  the  jet  where  P^b  becomes  the  surface  pressure  on  the  cone.  Conical  shock 
charts'^  indicate  that  for  M=2.9  and  0=25°,  the  surface  pressure  coefficient 

^.5  £l‘  =0.42  (3.1) 


where  p^  is  the  surface  pressure  of  the  cone  (therefore  P^b)’  Psi'^  freestream  static 
pressure,  and  qi  is  the  freestream  dynamic  pressure.  Thus, 


( 


0.42  ^ 

Po  P,i 


P,i- 


(3.2) 


To  obtain  the  3.25  pressure  ratio,  the  static  pressure  of  the  jet  must  be  3.25  times  P^b-  The 
pressure  ratio  (PR)  becomes  a  function  of  tunnel  stagnation  pressure  and  injector  stagnation 
pressure  given  by 


PR  =  1.68^. 
Pn 


(3.3) 


For  the  present  study,  p,j  =  386  kPa  and  Pt,  —  200  kPa. 


3.3  Data  Acquisition  Locations 

Both  mean  flow  and  hot-wire  data  was  acquired  at  two  downstream  locations  in  the 
wind  tunnel:  20  and  40  jet  exit  diameters  (x/d)  downstream  of  the  upstream  edge  of  the 
injector  nozzle. 

3.3.1  Mean  Flow  Data 

Mean  flow  data  was  taken  at  1 1  cross-stream  positions:  z/d  =  -5...0...+5.  At  each  z/d 
position  (referred  to  from  here  on  as  stations),  two  vertical  runs  had  to  be  made  because  of 
limited  run-time,  traverse  speed,  and  speed  limitations  for  accurate  probe  response.  Each 
run  covered  2.54  cm  of  the  flow  field,  from  y  =  -0. 16  to  -2.70  cm  on  the  first  run  and  from 
y  =  -2.70  to  -5.24  cm  on  the  second.  This  covered  approximately  13  y/d  stations.  See  Figure 


FIGURE  3.3:  Data  Acquisition  Areas 
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3.3.2  Hot-Wire  Data 

Single  overheat  hot-wire  data  was  taken  at  13  cross-stream  stations  from  z/d  =  - 
3...0...+3  at  every  half  station.  This  was  done  to  increase  the  resolution  in  the  region  around 
the  core  injector  flow.  A  single  vertical  profile  from  y  =  -0. 1 6  to  -3.97  cm  (10  y/d  stations) 
was  made.  Multiple  overheat  data  was  taken  only  at  the  center  (z/d  =  0)  downstream  station 
(x/d  =  40). 

3.4  Flow  Visualization 

Both  600  nsec  spark  schlieren  and  shadowgraph  photography  were  utilized  for  flow 
visualization. 

3.4.1  Shadowgraph  Setup 

Shadowgraphs  were  taken  of  the  entire  flowfield  encompassing  the  injector  and  both 
downstream  measuring  stations.  The  light  source  was  a  Cordin  Model  5401  arc  light, 
providing  a  600  nanosecond  spark.  Powering  the  spark  lamp  was  a  Cordin  Model  5205 
power  supply  set  at  approximately  5000  V  DC  for  both  schlieren  and  shadowgraph 
photography.  A  mirror  with  a  100  cm  focal  length  was  used  to  collimate  the  light  and  pass 
it  through  the  test  section.  All  photographs  were  taken  with  Type  52  Polaroid  film. 

3.4.2  Schlieren  Setup 

The  schlieren  setup  was  essentially  the  same  as  the  shadowgraph  setup  but  with  the 
addition  of  another  100  cm  focal  length  mirror  which  was  used  to  refocus  the  collimated 
beam  and  pass  it  through  the  knife  edge. 

3.5  Mean  Flow  Instrumentation 

The  wind  tunnel  was  manually  operated  by  first  opening  the  valve  to  the  vacuum 
tanks  then  the  valve  to  the  high  pressure  air  supply.  The  tunnel  was  shut  down  in  the  reverse 
order.  Tunnel  run  time  was  limited  primarily  by  the  available  vacuum  to  20-25  seconds. 


3.5.1  Settling  Chamber 

The  total  pressure  in  the  settling  chamber  was  held  constant  at  roughly  203  ±  4  kPa 
for  all  runs  and  was  recorded  on  the  data  acquisition  system  (DAS).  An  Endevco  0-690 
kPag  pressure  transducer  supplied  the  pressure  signal  which  was  also  used  as  the  DAS 
trigger.  A  mechanical  pressure  regulation  system,  consisting  of  a  Fairchild  pneumatic 
pressure  regulator  and  a  Leslie  pressure  reducing  valve,  kept  the  tunnel  pressure  constant 
during  the  run  (±4  kPa). 

As  indicated  previously,  total  temperature  was  monitored  by  an  Omega  Engineering 
thermocouple.  The  temperature  was  not  recorded  during  the  run  because  it  did  not  vary 
substantially  during  the  run  (298  ±3  K). 

3.5.2  Traverse  and  LVDT System 

The  traverse  was  comprised  of  an  Arrick  Robotics  MD-2  dual  stepper  motor  driver 
package  and  a  Size  23  Stepper  Motor.  The  stepper  motor  was  used  to  move  the  various 
probes  through  a  vertical  station  in  the  flow.  Translation  across  the  flow  was  performed 
manually  by  a  mechanical  screw/slide  assembly  holding  the  stepper  motor.  Maximum  speed 
of  the  stepper  motor  was  0.254  cm/sec.  The  position  of  the  traverse  was  monitored  by  a 
TransTek  Inc.  Model  0217  linear  voltage  displacement  transducer  (LVDT)  and  recorded  by 
the  DAS. 

3.5.3  Data  Acquisition  System 

Mean  flow  data  (tunnel  stagnation  pressure,  traverse  position,  and  probe  pressure) 
were  recorded  by  a  Nicolet  MultiPro  120  Digitizer^®  at  a  200  Hz  sampling  rate. 

3.5.4  Filters  and  Signal  Conditioners 

All  mean  flow  signals  were  conditioned,  filtered,  and  amplified  prior  to  acquisition 
by  the  DAS.  Signals  from  pressure  transducers  were  fed  to  Endevco  Model  4423  Signal 
Conditioners  which  supplied  excitation  voltage  to  the  transducers,  filtered  the  signal,  and 
amplified  it.  LVDT  signals  were  passed  through  a  TransTek  Model  1000-0012 
oscillator/demodulator. 
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3.5.5  Downstream  Pressure  Probes 

Mean  flow  parameters  were  measured  using  a  conventional  Pitot  probe  and  a  10° 
half-angle  cone-static  probe.  The  Pitot  probe  was  constructed  of  a  1.59  mm  outer  diameter 
stainless  steel  tube  inserted  into  a  3.175  mm  stainless  steel  tube  for  reinforcement.  An 
Endevco  0-103  kPag  transducer  recorded  both  Pitot  pressure  and  cone-static  pressure. 

The  cone-static  probe  consisted  of  a  10°  semivertex  angle  axisymmetric  cone.  The 
cone  tip  was  machined  from  stainless  steel  with  four  0.34  mm  (#80  drill)  pressure  taps  at 
90°  intervals  around  the  circumference  and  4.3  mm  from  the  tip.  The  four  pressure  taps  met 
in  a  common  chamber  where  they  were  averaged  to  account  for  misalignment  errors  of  up 
to  7°.“'  The  cone  was  soldered  to  1:59  mm  tubing  inserted  into  a  3.175  mm  sleeve  similar 
to  the  Pitot  probe.  The  pressure  taps  in  the  cone-static  probe  and  the  opening  in  the  Pitot 
probe  were  equidistant  from  the  probe  support  centerline. 

3.6  Hot-Wire  Instrumentation 

Data  acquisition  for  the  hot-wire  data  was  almost  identical  to  that  of  the  mean  flow 
data  with  the  substitution  of  two  hot-wire  channels  for  the  downstream  pressure  probe. 

3.6.1  Nicolet  Configuration 

A  total  of  four  channels  of  data  were  acquired  in  the  hot-wire  runs:  total  pressure, 
LVDT,  and  two  hot-wire  channels.  Total  pressure  and  hot-wire  channels  were  acquired  at 
10  kHz.  10  kHz  was  deemed  sufficient  to  capture  the  energy  of  the  turbulent  flow.  The 
probe  wires  were  tuned  to  provide  100-150  kHz  frequency  response.^^ 

3.6.2  Cross-Wire  Probes 

For  the  present  work,  TSI  Inc.  Model  1243-20  platinum  hot  film  probes  were  used. 
Each  wire  had  a  length  of  1  mm  and  a  diameter  of  5 1  pm,  with  a  coefficient  of  resistance 
of  0.24%/°  C  at  20  °C.  The  coefficient  was  assumed  constant  for  the  range  of  temperature 
encountered.  Model  1243-20  probes  were  used  to  acquire  “uv”  data  and  Model  1243 AN-20 
probes  were  used  to  acquire  “uw”  data.  The  probes  were  identical  except  for  their  wire 
orientation. 
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3.6.3  Anemometer 

A  TSI  Inc.  IFA  100  System  Intelligent  Flow  Analyzer,  constant  temperature 
anemometer  was  employed  for  the  present  study.  The  IFA  100  contained  a  Model  150 
anemometer  and  a  Model  1 57  Signal  Conditioner,  and  transducers  for  each  hot-wire  channel 
all  housed  in  the  Model  158  master  cabinet.  The  fFA  100  symmetrical  bridge  had  a  top 
resistance  of  50  Q  (RJ  and  allowed  for  an  external  bridge  arm  which  was  used  to  set  the 
overheat  ratio.  All  resistance  measurements  were  taken  with  a  Hewlett-Packard  3468A 
digital  multimeter  using  the  four-wire  method.  For  single  overheat  data,  the  external 
resistance  was  set  so  the  overheat  ratio  was  2.03±0.04.  For  multiple  overheat  data,  the 
external  resistance  was  varied  based  on  the  resistance  the  individual  wire  according  to  the 
following: 


Overheat  Number 
1 

2 

3 

4 

5 

6 

7 

8 


Overheat  Ratio 
1.25 
1.35 
1.50 
1.60 
1.75 
1.85 
1.95 
2.10 


For  each  overheat  ratio,  a  data  run  and  a  calibration  run  was  made.  While  some  fidelity  in 
matching  pressures  from  run  to  run  is  lost,  great  gains  are  made  in  spatial  resolution  over  the 
overheat  scanning  circuit  method  of  Bowersox. 


3. 7  Computational  Facilities 

Data  acquisition  was  accomplished  with  both  a  Zenith  486  and  a  Compaq  386  PC. 
Traverse  operation  was  performed  by  a  Zenith  286  and  a  Compaq  386  PC  (computers  were 
upgraded  during  the  test  sequence).  Data  reduction  was  performed  on  486  PCS  and  AFIT’s 
Sun  Sparc  20  work  stations.  Two-dimensional  plots  were  generated  with  Grapher  for 
Windows  and  contour  plots  were  generated  with  Mathcad  5.0+. 
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IV.  Data  Reduction  Techniques 


This  chapter  describes  the  data  reduction  techniques  for  both  the  mean  flow  and  hot¬ 
wire  (both  single  and  multiple  overheat)  data. 

4. 1  Mean  Flow 

All  mean  flow  data  were  combined  to  obtain  flowfield  Mach  number.  For  each  run, 
total  temperature  (Tn),  tunnel  stagnation  pressure  (pt,),  and  Pitot  pressure  (p,2)  or  cone-static 
pressure  (p^)  was  recorded.  The  Mach  number  can  be  calculated  by  the  ratio  of  cone-static 
to  Pitot  pressure  in  the  ffeestream.  This  ratio  as  a  function  of  Mach  number  for  the  10°  cone 
was  determine  by  a  least-square  fit  by  Bowersox;® 

—  =  -0.052976  +  4.6840  X  -  1 8.6786  +  50.7006  x"'  -54.1577x^  (4.1) 

M 


where  x  =  pc/pa-  Equation  (4.1)  is  valid  for  Mach  numbers  between  1.5  and  4.4  with  a 
standard  deviation  of  0.06%. 

4.2  Turbulent  Flow 

The  data  reduction  methods  for  the  fluctuating  component  of  the  flowfield  is 
considerably  more  difficult  than  the  mean  flow  component.  This  section  outlines  the  data 
reduction  methods  used  for  both  single  and  multiple  overheat  cross-wire  anemometry. 

4.2.1  General  Theory 

For  turbulent  compressible  flow,  the  Nusselt  number  of  a  cylinder  (ie  hot-film  probe) 
has  the  following  functional  relationship^^ 

Nu  =  Nu{Lld,  M,  Pr,  Re^,  t)  (4-2) 


where  L/d  is  the  wire  aspect  ration,  M  is  the  Mach  number,  Pr  is  the  Prandtl  number,  Re^  is 
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the  effective  Reynolds  number  based  on  wire  diameter,  and  t  is  the  temperature  loading 
factor  [t  =  “  T^)  /  is  the  wire  temperature  and  is  the  equilibrium  temperature 

the  wire  would  obtain  if  placed  unheated  in  the  flow.  Eqn.  (4.2)  reduces  to 

Nu  =  Nu(Re^,  t)  .  (4.3) 


for  the  following  conditions:  high  Reynolds  number  (Re  >  20),  constant  Prandtl  number 
flows  where  the  Mach  number  normal  to  the  wire  is  greater  than  approximately  1.2,  and 
where  the  wire  aspect  ratio  is  large  (L/d  »  1).“  For  wires  normal  to  the  flow,  hot-wire  data 
has  been  found  to  collapse  onto  the  following  curve 

Nu  -  a  sj^  +  b  (4-4) 


When  using  (4.4),  the  hot-wire  must  be  calibrated  at  each  wire  temperature. 

The  Nusselt  number  is  also  proportional  to  the  power  supplied  to  the  wire  (q  J 


Nu  = 


<1. 


(4.5) 


where  and  i,,,  =  V,  /  (R.,  +  R,  +  Rl)  (iw=wire  current,  R^=wire  resistance,  V,=wire 

voltage).  The  Nusselt  number  can  then  be  expressed  as  the  following,  if  T^  is  assumed  to 
be  T,  in  both  the  calibration  and  data  reduction  process 


Nu 


H'  W 


(4.6) 


To  arrive  at  the  appropriate  mean  and  fluctuating  equations,  the  voltage,  Reynolds 
number,  and  total  temperature  are  replaced  by  their  mean  and  fluctuating  components,  apply 
the  binomial  theorem,  and  noting 


n 


(4.7) 
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where 


(R  +R  +R,f 

r  =  Tii  k 

^0  ^  ^0 


(4.8) 


and  n,;  =  0.89.  The  detailed  procedure  to  derive  (4.7)  is  presented  in  Bowersox.^'^  Solving 

for 


V  ' 

{  Reo  '] 

(  'f’/'l 
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— 

V 

Re  o 

T, 

w 

\  ) 

\  ‘ 

(4.9) 


where  the  hot-wire  sensitivities  (f  andg)  are  given  as  the  following 


and 


(4.10) 


From  (4. 10),  it  is  clear  that  ^Re^  and  must  be  known  in  order  to  evaluate  the  hot-wire 
sensitivities.  Equation  (4.7)  can  be  rewritten  as 

(4.11) 


where  .r,  =  TJ,  and  -  =  vlj(C,a.  TJ-kJa..  Within  (4. 11),  /  is 

the  overheat  number  and  C  is  Cq  where  ko  is  replaced  with  k^  As  a  minimum,  two  overheat 
ratios  are  necessary  to  determine  the  two  unknowns  and  .  If  more  than  three 
overheat  ratios  are  used,  then  a  least  squares  analysis  can  be  used  to  obtain  ^Rs^  and  T^. 

Squaring  then  averaging  Eqn.  (4.9)  yields  the  following  hot-wire  fluctuation  equation 


fi 


^  Re  o  '  ^  ^ 


Re  o 


^f,S. 


^  Re  o'  T'^ 

e  i 


.  Re  o  T,  , 


I  T'\ 


+  e 


s, 


f. 


(  ..  A 


(4.12) 


The  general  least  square  (GLS)  method  yields  the  following  3x3  system,  where  i  again 
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denotes  the  overheat  number  and  where  the  summations  are  over  N,  the  number  of  overheat 
ratios  used. 
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It  is  important  to  note  that,  for  the  GLS  method,  any  errors  introduced  in  the  experiment 
through  scatter  are  amplified  by  a  combined  power  of  four.  This  inherent  disadvantage  of 
the  GLS  approach  led  to  the  development  of  a  new  method,  the  quadratic  least  square  (QLS) 
method.^ 

The  QLS  method  begins  by  assuming  that  the  error  in  the  measured  hot-wire  voltages 
are  purely  random  and  that  all  bias  errors  are  contained  in /and  Defining  t*  as 


T*  = 


(4.14) 


and  rewriting  the  hot-wire  sensitivities  as 
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(4.12)  can  be  rewritten  as 
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A  minimum  or  three  overheat  ratios  are  necessary  to  solve  for  the  unknowns  in  (4.16).  If 
more  than  three  are  used,  then  the  QLS  method  can  be  used  to  evaluate  the  right-hand  side 
of  (4. 1 6).  The  turbulence  results  can  then  be  obtained  by  solving  (4. 1 7). 

It  is  important  to  note  that  any  error  due  to  experimental  scatter  is  amplified  by  a 
combined  power  of  two.  The  GLS  method  and  the  QLS  method  yield  identical  results  if 
(4. 1 5)  is  valid  and  if  the  experimental  scatter  in  the  hot-wire  sensitivities  is  purely  random; 
the  QLS  method  only  reduces  the  effects  of  bias  error  in  f  and 
The  next  step  is  to  specify  Re^  for  the  cross-wire  probe. 


4.2.2  Cross- Wire 

To  obtain  cross- wire  turbulence  results,  Re^  must  be  related  to  the  x  and  y  tunnel 
coordinates.  The  following  analysis  is  presented  in  full  in  Bowersox.^^ 

It  can  be  shown  that 

Re  ’  =  Re^  +  "  Re^  =  A ,  Re"^  +2A^  Re^ Re^  +  A^  Re^  (4. 1 8) 


where  Aj  are  given  by  the  following 
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(4.19) 


=  cos“((I))  +  kl  sin‘((()) 
A^={\  -kl)  cos((}))  sin((j)) 
=  k^  cos^(4))  +  sin“((j)) 


and  (j)  is  the  hot-wire  incidence  angle. 

Replacing  Re^,  Re^,  and  Re^.  by  their  mean  and  fluctuating  components,  it  can  be 
shown  that 
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where  B,  =  A1/B3,  B2  =  A2/B3,  B,  =  A,  +  2A2R<,,  Ro  =  p  v/p  ?/,  and,/  is  the  wire  number  on 
the  cross-wire  probe.  Equation  (4.20)  assumes  Rq'  «  1 . 


4.2.3  Single  Overheat  Ratio  Technique 

For  many  compressible  flows,  the  total  temperature  fluctuations  are  small;  T,  is 
essentially  constant  and  d(T^  =  0  (adiabatic  flow).  For  adiabatic  flow  with  Pr  =  1,  provided 
T,'  =  0,  the  following  results 
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But  for  air,  Pr^^l  and  Tt’^0  even  for  flows  that  are  adiabatic  in  the  mean.  Therefore,  a 
correction  factor  can  be  added  to  account  for  Pr^  1  and  Tj’  ^^0 
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Typically,  k  e  [0.3,  0.5],  but  k  =  1 .0  was  used  in  this  study  based  on  the  high  overheat  ratio 
used  (T,’-0).  Defining  0s-k(y-1)M^  and  assuming  Ro^0  «  1-0,  then  it  can  be  shown  that 
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The  total  temperature  fluctuation  equation  (4.27)  now  becomes 
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where  a  =  1.0  /  [1  +0.5(y  -  l)M^]  and  P  =  (y  -  1)“^^  ■  With  the  above  equations  (4.22)- 
(4.25),  the  single  overheat  cross-wire  response  can  be  approximated  as 


'  Reo'^^ 


Re  o 


Reo^'  Reo^' 


Reo_  Reo^ 


(  13 

Re  o 

y 

2 

f  V 

(  V  '1 

w 

Re  o 

\  y 

V 

\  »' 

1 

V 

\  ^ ) 

2; 

V 

\  «’/ 


(  V 


V 

V  »'/ 


(4.26) 


where  S  is  the  single  overheat  matrix  whose  coefficients  are  given  below  (note  that  the 
coefficients  become  greatly  simplified  with  the  assumption  Rq-O) 
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Past  work  has  shown  for  boundary  layers  RyxE  [0.6,  0.9].  There  is  little  data,  however, 
concerning  R^j,  Rvt,  and  R,,,^  for  free  mixing  layers.  The  present  study  invoked  the  typical 
assumption  that  T,’=:0  (the  second  term  on  the  right-hand  side  of  Eqn.  (4.9)  was  neglected). 
This  arises  from  setting  k  =  1  or  R„,t  =  0  and  is  generally  accepted  to  be  accurate  for 
overheat  ratios  greater  than  1.9  (all  single  overheat  data  in  this  were  taken  at  OHR~2.0). 
This  is  true  because  the  sensitivity  g  decreases  as  the  OHR  is  increased.  Thus,  values  for 
R^.  1 ,  R,  I ,  and  R,,,t  were  not  needed  here. 
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All  of  the  single  and  multiple  overheat  ratio  data  reduction  techniques  have  been 
coded  in  the  FORTRAN  program  MSHeAR.^^ 

4.2.4  Separation  of  Turbulence  Variables 

The  information  provided  by  multiple  overheat  cross-wire  anemometry  in  supersonic 
flow  includes  the  following  “conservative”  variable  turbulence  data 

(P «,)'  (P ")'  (p »,)'  t;  (4.29) 


In  order  to  separate  the  “conserved”  turbulence  data  into  “nonconservative”  variables,  two 
assumptions  were  made. 

The  first  assumption,  that  the  pressure  fluctuations  are  small  compared  to  the  density 
and  temperature  fluctuations,  is  still  controversial.  Kistler^'’  has  suggested  that  p’  is 
proportional  to  u’^,  which  is  second  order  and  neglected.  The  validity  of  the  p’~  0(2)  (or 
p’  ~  0)  assumption  has  been  experimentally  verified  for  a  Mach  4.0  free  mixing  layer.  '^The 
second  assumption  is  the  thermally  perfect  gas  assumption  (p  =  pRT).  With  first-order 
assumptions,  the  following  “separated”  results  can  be  obtained,'^  which  include  p’ 
information  for  completeness 
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where  a  =  1.0  /  [1  +0.5(y  -  1)M^  and  P  =  (y  -  l)ccM^.  With  Equations  (4.30)-(4.32)  and  the 
p’  =  0  assumption,  all  the  turbulent  shear  terms  in  the  RANS  can  be  obtained. 
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V.  Results 


This  chapter  presents  the  results  of  the  low  angle  injection  compressible  turbulence 
data  acquisition  and  reduction.  Nominal  tunnel  parameters  (without  injection)  were 
determined  to  be  =  2.9,  P(„  =203  ±  4  kPa,  and  T,„  =  298  ±3  K  (see  Table  5. 1  for  a  full 
summary  of  mean  flow  parameters).  Shadowgraphs  of  the  entire  flowfield  and  schlierens 
of  the  area  of  injection  were  taken.  Mean  flow  and  single  overheat  data  were  taken  at 
numerous  stations  across  the  flow  at  two  locations  downstream  of  the  injector.  Multiple 
overheat  data  was  taken  at  the  furthest  downstream  station  (x/d  =  40)  at  the  centerline  (z/d 
=  0). 


Table  5. 1  Summary  of  Mean  Flow  Parameters 


M„ 

2.9 

P.CO 

200  kPa 

T,. 

298  K 

Re/m„ 

15  X  10" 

1.8 

P.1 

386  kPa 

T., 

298  K 

Re/mj 

25  X  10" 

Uj/U- 

0.86 

P/P:. 

3.29 

(pu)/(pu)„ 

2.82 

P/P» 

5.44 

1.65 

5.  /  Shadowgraph  and  Schlieren  Analysis 

Both  shadowgraphs  and  schlierens  provide  useful  qualitative  data  on  the  flowfield. 
The  exposure  time  (0.6  psec)  was  too  long  to  freeze  the  eddies  which  would  have  traveled 
0.5  mm  during  the  exposure.  Nonetheless,  unusual  features  of  the  injected  flow  are  visible. 

Figure  5.1  is  a  shadowgraph  of  the  entire  test  section.  The  cone-static  probe  is 
visible  at  the  x/d  =  40  station.  The  strong  shock  produced  by  the  injection  is  clearly  visible. 
The  weak  crossing  shocks  in  the  center  of  the  test  section  below  the  injector  are  formed  by 
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slight  misalignment  of  the  tunnel  sections.  The  injector  shock  can  be  seen  reflected  from 
the  lower  wall  and  boundary  layer.  This  shock  and  its  reflection  are  clearly  visible  in  the 
mean  flow  Mach  number  contours  (Figure  5.3).  The  white  line  at  the  bottom  of  the 
shadowgraph  is  simply  a  reflection  from  the  light  source  off  the  tunnel  floor. 

Figure  5.2  is  a  schlieren  of  the  first  half  of  the  flow  field.  The  boundary  layers  are 
clearly  visible.  Shown  in  greater  detail  is  the  injection  point  and  associated  flow  structures. 
An  unusual  feature  of  this  picture  is  the  uneven  shock  downstream  of  the  injector  shock. 
This  shock  is  somewhat  unsteady  as  it  was  shaped  slightly  differently  in  various  schlierens 
and  shadowgraphs.  It  is  most  likely  a  combination  of  reattachment  shock  as  the  boundary 
layer,  initially  separated  by  the  jet,  reattaches  to  the  tunnel  ceiling  and  a  shock  from  the 
downstream  lip  of  the  injector. 

5.2  Conventional  Probes 

Conventional  probe  data  (Pitot  and  cone-static)  were  taken  at  both  x/d  stations.  The 
flow  area  mapped  out  is  approximately  5.08  x  3.96  cm  (nearly  50%  of  the  total  cross 
section)  (see  Fig.  3.3).  Normalized  Pitot  and  cone-static  pressures  as  well  as  Mach  number 
contours  are  provided  in  Figs.  5. 3-5. 5.  Approximately  40  points  of  raw  data  were  averaged 
per  point  plotted  (approximately  0.5  mm  grid  spacing).  Eleven  stations  across  the  flow  were 
taken  for  a  grid  spacing  of  about  4  mm.  In  addition  to  providing  insight  into  the  mean  flow 
character  of  the  flow,  the  mean  flow  data  had  two  additional  purposes;  to  provide  needed 
input  for  the  hot-wire  data  reduction  and  to  show  the  extent  of  jet  core  penetration  to  narrow 
the  data  taking  for  the  hot-wire  data. 

No  exact  definition  of  jet  core  size  exists,  but  one  can  determine  the  extent  to  which 
the  jet  penetrates  the  freestream  by  measuring  the  position  of  the  jet  center.  One  can 
approximate  the  outer  dimensions  of  the  jet  by  mapping  Mach  number,  shear  stress, 
turbulence,  etc. 

5.3  Conventional  Probe  Results 

5.3.1  Station  x/d=20 

Based  on  the  Mach  contour  (Fig.  5.3)  the  jet  core  penetrated  to  y/d  =  -3.0.  The  jet 
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FIGURE  5.2;  Schlieren  of  Injector  Flowfield 
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is  roughly  symmetrical  and  encompasses  z/d  e  [-1.6, 1.6],  The  average  diameter  of  the  jet 
is  nominally  3. 2d.  The  boundary  layer  on  the  floor  of  the  tunnel  is  approximately  1.5d, 
while  the  boundary  layer  on  the  ceiling  (injector  side)  of  the  tunnel  has  increased  to  about 
2.0-2. 5d  because  of  the  action  of  the  jet.  It  is  difficult  to  determine  the  exact  thickness  of 
the  boundary  layer  at  this  point  because  of  the  interaction  with  the  jet. 

5.3.2  Station  x/d=40 

Again,  based  on  the  Mach  contour  (Fig.  5.3)  the  jet  core  has  penetrated  to  y/d  =  -4.5 
with  an  average  diameter  of  roughly  3.6d.  The  jet  has  expanded  only  slightly  to  encompass 
z/d  6  [-1 .8,  1 .8].  The  boundary  layer  has  not  grown  significantly  by  this  point  because  of  the 
short  distance  between  stations  (about  7.9  cm).  This  can  be  seen  in  Fig.  5.2  and  to  a  lesser 
extent.  Fig.  5.1.  The  jet  exhibits  only  slight  asymmetry.  There  are  no  strong  shocks 
intersecting  the  upper  boundary  layer  between  the  injector  and  x/d  =  40.  The  reflection  of 
the  jet  shock  (although  not  as  sharply  defined)  is  at  y/d  =  -10.  Clearly  visible  at  z/d  =  3  (and 
to  a  much  lesser  extent  at  z/d  =  -3),  is  a  weak  vertical  structure  which  is  thought  to  be  a 
shock  reflected  off  the  sidewall  of  the  tunnel.  This  weak  shock  most  likely  arises  from  a 
slight  misalignment  of  tunnel  sections.  The  Mach  number  change  across  the  weak  shock  is 
less  than  0.05.  There  is  a  complex,  three-dimensional  shock  pattern  in  the  test  section. 
While  the  shock  visible  in  Figs.  5. 1  and  5.2  reintersects  the  tunnel  ceiling  well  downstream 
of  the  areas  of  interest,  the  same  cannot  be  said  for  the  sidewall  shocks  which  have  only  half 
the  distance  to  travel  from  the  injector.  Another  possibility  is  that  of  another  set  of  weak 
misalignment  shocks  from  the  nozzle/test  section  interface.  The  normalized  Pitot  and  cone- 
static  pressures  (Figs.  5.4  and  5.5)  follow  the  same  general  contour  as  the  Mach  number 
contour.  Based  on  the  jet  core  penetration  from  Fig.  5.3,  it  was  decided  to  take  hot-wire  data 
in  the  area  z/d  e  [-3,3]  and  y/d  e  [0,9.8]. 

5.4  Cross-Wire  Probes 

Single  overheat  cross-wire  data  were  obtained  at  every  half  z/d  station  from  z/'d  e  [- 
3,3]  and  at  x/d  =  20  and  40.  The  mean  flow  pressure  data  were  linearly  interpolated  between 
stations  as  needed  for  the  hot-wire  data  reduction.  Multiple  overheat  data  was  taken  only 
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at  the  x/d  =  40  centerline  station. 


5.4.1  Calibration  Technique 

The  cross-wire  probes  were  calibrated  over  the  range  of  Reynolds  numbers 
anticipated  in  the  flow  field,  approximately  15  -  25x  10^/m,  based  on  the  calculated  extremes 
in  the  tunnel  and  jet  (see  Table  5.1).  The  minimum  stagnation  pressure  at  which  the  tunnel 
could  operate  was  approximately  159  kPa  and  the  maximum  attainable  by  the  pressure 
regulation  and  supply  system  was  290  kPa.  The  anemometer  was  tuned  in  accordance  with 
the  IFA  1 00  manual  at  roughly  200  kPa  whenever  a  probe  was  changed.  All  calibrations 
were  performed  at  the  approximate  center  of  the  test  section  with  the  injector  off  For  the 
calibration  run,  the  tunnel  stagnation  pressure  was  varied  from  the  minimum  to  as  close  to 
the  maximum  as  was  possible.  Roughly  ten  data  points  were  taken  from  the  calibration 
process  to  produce  the  variables  a  and  b  in  Eqn.  (4.10). 

For  single  overheat  data,  a  calibration  run  was  made  at  the  beginning  and  end  of  a 
data  acquisition  sequence.  Generally,  all  thirteen  z/d  stations  were  taken  in  a  2-3  hour 
period.  The  calibration  runs  were  compared  to  determine  if  any  electronic  drift  or  probe 
contamination  had  occurred.  Neither  were  noted  during  any  test  sequence.  For  multiple 
overheat  data,  a  calibration  run  was  performed  for  every  overheat  ratio.  The  sequence 
involved  setting  the  overheat,  taking  the  calibration  data,  then  taking  the  flow  data.  This 
sequence  was  repeated  eight  times  per  probe  type. 

5.4.2  Single  Overheat  Data 

Two  probe  types  were  used:  the  TSl  1243-20  (“uv”  probe)  and  the  TSI  1243AN-20 
(“uw”  probe).  Probe  wire  resistances  as  measured  by  the  four-wire  method  and  overheat 
ratios  used  appear  in  Table  5.2.  All  single  overheat  ratios  were  kept  in  a  narrow  band  OHR 
6  [1.99,  2.06].  At  high  overheat  ratios  the  probes  are  less  sensitive  to  temperature 
fluctuations  (i.e.  the  second  term  of  the  right-hand  side  of  Eqn.  (4.9)  is  small  compared  to 
the  rest).  That  is,  the  wire  temperature  is  much  greater  than  the  flow  temperature. 
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Table  5.2  Cross-Wire  Probe  Operating  Resistances  and  Overheat  Ratios 


Probe 

Station 

Wire  1,  Q  (OHR) 

Wire  2,  Q  (OHR) 

uv 

x/d=20  (SOH) 

5.90  (2.00) 

5.94(1.99) 

uw 

x/d=20  (SOH) 

5.59  (2.02) 

6.34  (2.01) 

uv 

x/d=40  (SOH) 

5.95  (2.06) 

5.99  (2.04) 

uv 

x/d=40  (MOH) 

5.91  (variable) 

5.98  (variable) 

uw 

x/d=40  (SOH) 

6.28(2.05) 

6.03  (2.03) 

uw 

x/d=40  (MOH) 

6.20  (variable) 

5.96  (variable) 

A  full  set  of  data  was  not  taken  at  x/d  =  20  with  the  “uw”  probe.  A  series  of  probe 
wire  failures,  probably  the  result  of  faltering  filters  allowing  small  particulates  into  the  flow, 
left  only  stations  z/d  =  -3,  -2,  -1,  and  0  measured  (half  stations  were  not  measured). 
Therefore,  the  spatial  resolution  of  any  data  at  x/d  =  20  requiring  the  “uw”  probe  is  only  half 
as  good  as  the  other  stations.  All  contour  plots  using  the  “uw”  probe  at  x/d  =  20  have  the 
left  half  reflected  onto  the  right.  That  is,  station  z/d  =  3  was  used  as  z/d  =  -3,  etc.  This  is 
a  fair  practice  since  with  only  one  exception  (flow  angles.  Fig  5. 13  and  velocity  component 
w.  Fig.  5.50)  the  jet  exhibits  a  high  degree  of  symmetry. 

5.4.3  Multiple  Overheat  Data 

Eight  overheat  ratios  were  taken  with  each  probe  type  at  overheats  based  on  the  chart 
in  Section  3.6.3.  The  overheat  ratios  ranged  from  1 .25  to  2. 1,  but  the  highest  overheat  and 
lowest  ratios  were  generally  not  used  because  they  exhibited  instabilities.  Continuous 
multiple  overheat  data  were  taken  (one  overheat  ratio  per  run)  and  combined. 
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Table  5.3  Actual  Overheat  Ratios  for  Multiple  Overheat  Data 


“uv” 

Wire  1 

Probe 

Wire  2 

“uw” 

Wire  1 

Probe 

Wire  2 

Overheat 

1 

2.08 

2.06 

t 

t 

2 

1.99 

1.97 

1.96 

1.96 

j 

1.83 

1.80 

1.83 

1.82 

4 

1.73 

1.71 

1.75 

1.73 

5 

1.63 

1.61 

1.58 

1.61 

6 

1.54 

1.52 

1.51 

1.48 

7 

1.35 

1.33 

1.34 

1.36 

8 

1.26 

1.24 

1.25 

1.27 

^  The  probe  failed  on  this  run. 


5.5  Cross-Wire  Probe  Results 

The  “uv”  and  “uw”  cross-wire  probes  generally  showed  good  fidelity  on  quantities 
which  both  probes  measured.  Turbulence  intensity  (Figs  5.6  and  5.10)  was  within  ±2%  based 
on  the  root-mean-square  (rms)  quantities.  Favre  terms  (Figs  5.23  and  5.24)  matched  to 
within  ±0.001.  Separated  turbulence  intensities  (Figs.  5.29,  5.30,  5.31,  and  5.32)  were 
within  ±1%  rms.  Density-velocity  correlations  (Figs.  5.37  and  5.38)  were  within  ±0.002. 
Finally,  mean  velocities  (Figs.  5.47  and  5.48)  were  within  ±10  m/sec.  One  exception  to  the 
previous  generality  is  in  measurements  in  the  boundary  layer  (upper  reaches  of  Figs.  5.53(a) 
and  5.54(a)). 

This  difference  is  thought  to  arise  out  of  both  the  scale  and  orientation  of  the  probes. 
Both  probes  are  easily  larger  than  the  smaller  eddies  encountered  in  the  flow  field  and  would 
be  expected  to  evince  some  scatter  between  them.  More  importantly  perhaps,  is  the 
orientation  of  the  probe  wires.  The  “uv”  probe’s  wires  are  oriented  vertically 
(perpendicular  to  the  tunnel  floor).  In  this  orientation,  the  probe  wires  are  less  affected  by 
sharp  vertical  gradients  as  are  found  in  the  boundary  layer.  The  “uw”  probe’s  wires  are 
oriented  horizontally  so  that  in  a  sharp  vertical  gradient  the  two  wires  could  experience 
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vastly  different  flows.  One  solution  to  the  problem  would  be  to  ensure  the  probe  wires  are 
never  perpendicular  to  such  a  gradient. 

5.5.1  Single  Overheat  Cross- Wire  Results 

The  single  overheat  results  are  shown  in  contour  plots  labeled  Figs.  5.6-5.52.  Each 
figure  contains  two  contour  plots  one  for  the  x/d  =  20  station  (a)  and  one  for  the  x/d  =  40  (b) 
station  except  for  Fig.  5.51  which  contains  a  single  plot  at  the  x/d  =  40  station. 

Perhaps  the  most  interesting  gross  feature  of  the  injected  flow  is  a  slight  bifurcation 
of  the  jet.  The  bifurcation  is  clearly  visible  in  Fig.  5.6  at  both  x/d  stations.  The  phenomenon 
is  also  apparent  in  all  turbulence  intensity,  Favre,  turbulent  kinetic  energy,  vorticity,  and 
most  shear  stress  graphs.  This  bifurcation  probably  arises  from  vortical  motion  in  the  jet 
itself,  not  unlike  the  wingtip  vortices  seen  on  aircraft.  No  similar  phenomenon  was  observed 
in  the  literature.  It  should  be  noted  that  most  literature  on  low-angle  injection  involves 
helium  injection  with  measurement  of  concentration  in  the  jet.  The  measurement  techniques 
employed  are  not  designed  to  measure  the  observed  vortical  motion. 

The  apparent  existence  of  the  bifurcation  is  bolstered  by  direct  flow  angle  (Figs.  5.9 
and  5. 13)  and  mean  flow  velocity  (Figs.  5.49  and  5.50)  measurements.  The  V  and  W  mean 
flow  velocities  were  combined  into  a  single  graph  of  vorticity  about  the  x-axis  (Fig.  5.51 ). 
The  figure  shows  the  bifurcations  to  be  two  counterrotating  vortices  of  nearly  equal  strength. 

Two  hypothesis  are  proposed  to  account  for  the  phenomenon.  First,  the  bifurcation 
may  arise,  or  at  least  be  bolstered,  by  the  design  of  the  injector.  Recall  from  Fig.  3.2  and 
Section  3.2.2  that  the  lower  half  of  the  injector  settling  chamber  is  spherical.  This  could 
induce  opposing  vortices.  However,  boundary  layer  vorticity  is  the  more  likely  player. 

As  the  jet  enters  the  boundary  layer  it  causes  a  temporary  separation  which  serves 
to  entrain  the  boundary  layer  into  the  jet  (see  Fig.  5.3).  Additionally,  flow  around  the  jet 
would  lead  to  some  vortical  motion.  The  jet  was  modeled  as  a  half  cone  in  the  flow  for  the 
purpose  of  determining  which  is  true  initially,  but  as  the  jet  rises  out  of  the  boundary 
layer  it  behaves  more  like  a  cylinder  inclined  to  the  flow.  As  the  flow  wraps  around  the 
cylinder,  it  can  impart  some  vortical  motion  on  the  jet  itself 

In  Fig.  5.6,  and  subsequent  figures,  the  jet  appears  more  rounded  than  in  the  Mach 
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contour.  While  the  measured  center  of  penetration  of  the  core  does  not  change  appreciably 
(<  0.  Id),  the  extent  of  influence  of  the  jet,  which  one  could  consider  the  diameter  of  the  Jet, 
has  increased.  At  x/d  =  20,  the  Jet  diameter  is  4.0d  and  at  x/d  =  40,  it  is  6.0d.  Measures 
taken  from  other  turbulence  intensity  plots  and  the  turbulent  kinetic  energy  plot  (Fig.  5.14) 
show  roughly  (±0.2d)  the  same  Jet  diameter.  Turblulence  intensity  is  likely  a  better  indicator 
of  Jet  size  in  the  freestream. 

As  shown  in  previous  studies,'^  compressibility  accounts  for  about  50-75%  of 
Reynolds  shear  stress.  Figs  5. 15-5. 18  show  the  separated  Reynolds  shear  components  as 
measured  by  the  “uv”  probe.  Fig.  5. 15  shows  the  incompressible  portion  (first  term  of  Eqn. 
(2.14)),  while  Figs.  5.16  and  5.17  show  the  compressible  components  (second  and  third 
terms  of  Eqn.  (2. 14)).  Note  the  maximum  values  of  Figs.  5. 15-5. 17.  The  incompressible 
portion  is  approximately  half  that  of  the  first  compressible  term  while  the  second 
incompressible  term  is  an  order  of  magnitude  smaller  than  the  first.  The  third  term  can  be 
safely  ignored  in  turbulence  models,  however,  since  the  second  term  accounts  tor  roughly 
67%  of  the  total  Reynolds  shear  stress  (Fig.  5. 18)  it  cannot  be  ignored  and  must,  in  fact,  be 
rigorously  included.  The  same  analysis,  with  the  same  result,  can  be  applied  to  the  Reynolds 
shear  measured  by  the  “uw”  probe  (Figs.  5.19-5.22). 

Since  most  modem  turbulence  models  are  based  on  the  Favre-averaged  Navier- 
Stokes  equations  Favre  data  is  included  as  Figs.  5.23-5.28.  Notice  that  the  Favre  shear  (Figs. 
5.27  and  5.28)  match  nearly  exactly  the  incompressible  portion  of  the  Reynolds  shear  (Figs. 
5. 1 5  and  5. 19).  This  result  is  expected  since  the  FANS  equations  bear  a  close  resemblance 
to  the  incompressible  RANS. 

For  completeness,  the  separated  turbulence  intensities  (pVp)™^,  (u7u)n„„  (v7u)rn,„ 
and  (w7u)™,  are  included  as  Figs.  5.29-5.34.  Also,  the  velocity-velocity  and  velocity- 
density  correlations  are  included  as  Figs.  5.35-5.40. 

Turbulent  shear  stress  obtained  through  the  turbulence  transformation  process 
discussed  in  Section  2.4  are  presented  as  Figs.  5.41-5.46.  For  data  obtained  with  the  “uv” 
probe.  Fig.  5.41  shows  the  first  term  of  Eqn.  (2.21),  Fig.  5.42  shows  the  second  term,  and 
Fig.  5.43  shows  the  total.  Figs.  5.44-5.46  show  the  “uw”  probe  data.  The  data  indicate  the 
first  term  comprises  approximately  90%  of  the  total  and,  subsequently,  the  second  term  is 
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an  order  of  magnitude  smaller  than  the  first.  This  result  is  important  because  it  shows  the 
first  term  to  be  of  primary  importance  for  modeling  purposes.  It  also  indicates  the  validity 
of  the  p’=0  assumption  discussed  in  Section  4.2.4  does  not  have  a  significant  impact  on  the 
total  transformed  shear.  This  arises  because  only  the  second  term  of  Eqn.  (2.21)  is  reliant 
on  separation  of  variables  which  requires  the  p’=0  assumption.  Thus,  even  if  the  p’=0 
assumption  is  found  to  invalid,  the  total  transformed  shear  would  be  minimally  affected. 

Mean  velocities  and  mean  mass  flux  are  provided  in  Figs.  5.47-5.50  and  Fig.  5.52 
respectively.  The  data  from  Figs.  5.49  and  5.50  were  used  to  generate  vorticity  plot  (Fig. 
5.51). 

5.5.2  Multiple  Overheat  Cross- Wire  Results 

Figs.  5.53-5.62  contain  the  multiple  overheat  data  taken  at  the  x/d  =  40  centerline 
station.  Figs.  5.53-5.56  show  comparisons  of  multiple  and  single  overheat  data  taken  at  that 
station. 

The  multiple  overheat  data  show  good  agreement  with  single  overheat  data  taken  at 
the  x/d  =  40  centerline  station.  Generally,  the  multiple  overheat  data  evinces  greater  scatter 
than  the  single  overheat.  This  is  true  because  the  multiple  overheat  data  was  collected  over 
as  many  as  eight  runs  while  the  single  overheat  was  collected  in  a  single  run.  Some 
difference  in  tunnel  conditions  from  run  to  run  is  inevitable.  Any  inherent  unsteadiness  in 
the  flow  will  exacerbate  the  scatter. 

Overall,  multiple  overheat  data  showed  “uw”  turbulence  intensities  (Fig.  5.54)  to  be 
less  than  that  measured  by  the  single  overheat  method  but  showed  “uv”  turbulence 
intensities  (Fig.  5.53)  to  be  greater.  The  opposite  trend  was  visible  in  shear  data  (Fig.  5.55). 
These  differences  are  within  the  expected  error  (see  Section  5.6). 

Flow  angle  data  (Fig  5.56)  showed  multiple  overheat  “uv”  data  to  be  in  close 
agreement,  however,  “uw”  flow  angle  data  show  the  single  overheat  flow  data  to  be  off  by 
-4.6°.  This  error  was  accounted  for  in  Fig.  5.13  by  adding  4.6°  to  the  entire  flowfield. 
Similar  analysis  of  the  r  component  of  the  mean  flow  indicated  an  error  of  approximately 
-50  m/sec.  Figure  5.50  was  corrected  as  well.  The  most  likely  source  of  this  error  is  a 
misalignment  of  the  “uw”  probe.  The  nature  of  hot-wire  anemometry  tends  to  reduce  small 
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angular  errors,  but  only  when  these  errors  are  consistent  in  both  the  calibration  and  data 
acquisition.  Changing  the  angle  of  the  probe  to  the  flow  ±5°  in  the  data  reduction  only 
varied  the  final  determination  of  flow  angle  by  ±0.5°.  However,  calibrating  the  probe  at  one 
angle  and  running  the  data  reduction  at  a  different  angle  produces  a  change  in  flow  angle 
nearly  equal  to  the  initial  difference.  It  is  very  likely  the  data  and  calibration  runs  were 
taken  with  different  probe  alignments. 

Total  temperature  turbulence  is  plotted  in  Fig.  5.57.  It  indicates  the  total  temperature 
fluctuated  a  maximum  of  approximately  3%  which  matches  the  maximum  difference 
between  single  and  multiple  overheat  data  in  Figs.  5.53-5.55.  The  assumption  of  small  total 
temperature  fluctuation  appears  valid. 

Most  of  the  data  in  Figs.  5.58-5.62  reach  extrema,  maxima  or  minima,  near  y/d  =  - 
5.0.  This  may  indicate  the  center  of  the  jet  and  therefore  jet  penetration.  However,  as  seen 
in  the  single  overheat  data,  the  extreme  values  of  the  data  did  not  usually  occur  at  the 
centerline  station.  Usually,  the  extrema  occurred  at  the  center  of  the  bifurcations.  This 
leads  to  some  ambiguity  in  the  determination  of  jet  penetration  based  on  compressible 
turbulence.  Since  the  single  and  multiple  overheat  data  vary  to  some  degree,  more  multiple 
overheat  data  across  the  flowfield  will  be  necessary. 

Fig.  5.58(a)  shows  the  mass  flux/total  temperature  correlations  to  track  well  together 
until  the  boundary  layer/jet  interaction  area  is  reached  (y/d  e  [0,  4.5]).  Fig.  5.58(b)  indicates 
the  r  component  of  the  specific  turbulent  heat  flux  to  be  much  less  than  the  other  two  except 
in  the  boundary  layer/jet  interaction  area.  This  makes  physical  sense  based  on  the  relative 
velocities  and  flow  angles.  At  the  centerline,  the  r  component  of  velocity  and  angle  would 
ideally  be  zero,  and  it  is  quite  close  to  the  ideal  over  most  of  the  centerline.  Figs.  5.13  and 
5.50  show  the  departure  of  the  flow  angle  and  velocity  from  the  ideal  nearer  the  wall. 

The  multiple  overheat  Reynolds  shear  stress  (Fig.  5.59)  complements  the  single 
overheat  data  showing  compressibility  to  account  for  60-70%  of  the  total  Reynolds  shear. 

Figs.  5.60-5.62  show  the  general  trends  discussed  thusfar;  extrema  at  y/d  =-5.0, 
increase  in  absolute  magnitude  of  "  components  in  the  boundary  layer/jet  interaction  area, 
and  general  departure  of  “uv”  and  “uw”  data  in  the  boundary  layer. 
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5.6  Error  Analysis 

Scatter,  or  error,  is  a  natural  part  of  taking  any  measurement,  in  particular  of  an 
unknown  quantity.  Presented  in  this  section  are  some  of  the  sources  of  error  in  this  work 
and  the  extent  to  which  these  errors  affect  the  accuracy  of  the  results.  The  L2  or  Euclidean 
norm"^  is  employed  as  the  measure  of  error  here.  The  L2  norm  is  given  by 


/  =  1 


where  /  indexes  the  various  errors  associated  with  measurement  of  the  quantity  x 
5.6.1  Conventional  Probes 

Volluz^'  indicates  turbulence  induces  about  ±690  Pa  error  for  both  Pitot  and  cone- 
static  probes.  The  calibration  error  is  approximately  500  Pa  for  both  probes.  The  settling 
chamber  pressure  has  an  uncertainty  of  about  1.7  kPa. 

The  probe  was  observed  to  flex  about  2.5  mm  while  the  tunnel  was  in  operation, 
therefore  the  probe  location  error  is  estimated  to  be  ±0.5%.  This  corresponds  to  a  Mach 
number  error  of  ±0.2  or  ±7%.  The  error  in  the  local  mass  flux,  based  on  previous  work,  is 
±2% 


5. 6. 2  Cross-  Wire  Probes 

Multiple  overheat  data  was  taken  continuously  and  averaged  over  every  2,048  points. 
Based  on  the  speed  and  data  acquisition  rate  of  the  probes,  the  position  error  is  0.5  mm, 
however,  the  largest  errors  are  associated  wnth  the  calibration.  Bowersox'^  indicates 
conventional  probe  errors  can  be  carried  over  into  the  hot-wire  analysis.  The  calibration 
errors  can  be  estimated  through  a  perturbation  analysis.  This  analysis  has  shown  error 
associated  with  the  probe  location  to  be  small,  primarily  because  the  probe  is  calibrated  at 
each  overheat  and  each  station.  As  shown  in  Section  5.5.2,  varying  the  angle  of  the  probe 
in  the  flow  as  much  as  5°  caused  only  a  0.5°  error  in  flow  angle.  Provided  the  probe  is  not 
removed  or  replaced,  misalignment  errors  will  be  minimized. 

The  error  analysis  of  Bowersox'^has  also  shovm  the  error  arising  from  linearization 


38 


of  the  hot-wire  response  equation,  through  the  binomial  theorem  (Chapter  4)  to  be  very 
small.  Bowersox*'  presents  the  following  hot-wire  errors  associated  with  similar  turbulence 
measurements  under  similar  flow  conditions  as  those  in  the  present  work; 

Table  5.4  Cross-Wire  Errors* 


e-r 

7.6% 

% 

3.6% 

Cu 

1.5% 

^pu 

2.1% 

e, 

1.0% 

eg 

1.0% 

^O'O’ 

8.4% 

*based  on  M=3.0 

Comparison  of  data  measured  by  both  the  “uv”  and  “uw”  probes  can  yield  insight 
into  the  errors  associated  with  the  present  flow.  As  stated  in  Section  5,5,  the  u  velocity  error 
is  10  m/sec  out  of  a  maximum  of  about  600  m/sec.  This  is  an  error  of  1 .5%  and  follows  well 
with  the  errors  shown  in  Table  5.4.  The  turbulence  intensity  error  is  roughly  0.02  rms. 
Maximum  error  in  the  Favre  data  is  approximately  12%.  The  separated  turbulence 
intensities  are  within  0.01  rms.  Density-velocity  correlations  have  a  maximum  12%  error. 

In  the  aggregate,  the  data  is  repeatable  within  12%  for  the  two  types  of  probes  which 
is  a  good  measure  of  the  repeatability  from  run  to  run  and  accounts  for  all  types  of 
measurement  errors:  probe  alignment,  temperature,  pressure,  voltages,  etc. 
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FIGURE  5.6:  Turbulence  Intensity  [(pu)7pu]^  (uv  probe) 
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FIGURES.?:  Turbulence  Intensity  [(pv)ypu]^  (uv probe) 
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FIGURE  5.8;  Shear  [(pu)'(pv)7(pu)^]  (uv  probe) 
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FIGURE  5.9:  Flow  Angle  [0pv/p„]  (uv  probe) 
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FIGURE  5.10:  Turbulence  Intensity  [(pu)7pu]^  (uw  probe) 
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FIGURE  5.11:  Turbulence  Intensity  [(pw)'/pu]^  (uw  probe) 
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FIGURE  5.15:  Separated  Reynolds  Shear  [-puV/pu^]  (uv  probe)  (xlO^) 
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FIGURE  5.16:  Separated  Reynolds  Shear  [-upV/pu^]  (uv  probe)  (xlO^) 


FIGURE  5.18:  Separated  Reynolds  Shear  [Reynolds  Shear]  (uv  probe)  (xlO^) 
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FIGURE  5.19:  Separated  Reynolds  Shear  [-pu'w'/pu^]  (uw  probe)  (xlO^) 
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FIGURE  5.20:  Separated  Reynolds  Shear  [-up'w’/pu^]  (uw  probe)  (xlO^) 


20.3 


59 


FIGURE  5.22:  Separated  Reynolds  Shear  [Reynolds  Shear]  (uw  probe)  (xlO^) 
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FIGURE  5.30;  Separated  Turbulence  Intensity  [p'/p]^  (uw  probe)  (xlO^) 
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FIGURE  5.31:  Separated  Turbulence  Intensity  [u'/u]^  (uv  probe)  (xlO^) 
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FIGURE  5.32:  Separated  Turbulence  Intensity  [u'/u]^  (uw  probe)  (xlO^) 
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FIGURE  5.33:  Separated  Turbulence  Intensity  [v'/uj^^  (uv  probe)  (xlO^) 
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FIGURE  5.34;  Separated  Turbulence  Intensity  [w'/uj^^  (uw  probe)  (xlO^) 
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FIGURE  5.35;  Velocity-Velocity  Correlation  [(u'v')/u^]  (uv  probe)  (xlO^) 
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FIGURE  5.36:  Velocity-Velocity  Correlation  [(u'w')/u^]  (uw  probe)  (xlO^) 


FIGURE  5.38:  Density- Velocity  Correlation  [(p'u')/pu]  (uw  probe)  (xlO^) 
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FIGURE  5.40:  Density- Velocity  Correlation  [(p'w')/pu]  (uw  probe)  (xlO^) 
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FIGURE  5.41:  Transformed  Shear  [(pu)'(pv)7(pu)^]  (uv  probe)  (xlO^) 
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FIGURE  5.42:  Transformed  Shear  [uv(p')7(pu)^]  (uv  probe)  (xlO^) 
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FIGURE  5.44:  Transformed  Shear  [(pu)'(pw)'/(pu)^]  (uw  probe)  (xlO^) 
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FIGURE  5.45:  Transformed  Shear  [uw(p')V(pu)^]  (uw  probe)  (xlO^) 
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FIGURE  5.49;  Mean  Velocity  [V  (m/s)]  (uv  probe) 
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FIGURE  5.51:  Vorticity  (coj  (rad/sec) 
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FIGURE  5.53:  Turbulence  Intensity  (uv  probe) 
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FIGURE  5,54:  Turbulence  Intensity  (uw  probe) 
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FIGURE  5.56:  Flow  Angle  (9) 
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FIGURE  5.57:  Turbulence  Intensity  [T  VT^  ] 
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FIGURE  5.58:  Mass  Flux/Total  Temperature  and  Specific  Turbulent  Heat  Flux 
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FIGURE  5.60:  Fa\  re  Variable  Data 
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FIGURE  5.62:  Velocity-Velocity  and  Density-Velocity  Correlations 


VI,  Conclusions  and  Recommendations 


6.1  Conclusions 

Mean  flow  data  and  flow  visualization  provided  insight  into  mean  flow  structure  and 
guidance  to  reduce  the  extent  to  which  the  flowfield  was  measured  with  hot-wire  probes. 

Single  overheat  anemometry  across  the  flowfield  at  each  station  revealed  turbulence 
intensity  and  turbulent  kinetic  energy  to  be  better  estimates  of  jet  core  penetration  and  jet 
size  than  mean  flow  data  alone.  The  single  overheat  data  also  confirmed  the  role 
compressibility  plays  in  the  Reynolds  shear  in  boundary  layer  and  free  mixing  layer  flows. 
This  study  determined  compressibility  to  account  for  approximately  67%  of  the  Reynolds 
shear  stress. 

In  particular,  the  data  indicate  the  last  two  terms  of  Eqn.  (2.14),  repeated  as 
Eqn.(6.t),  to  be  insignificant.  The  second  term  is  very  significant,  however,  and  is 
frequently  neglected  in  turbulence  models.  It  is  the  second  term  which  comprises  virtually 
all  the  compressible  portion  of  the  total  Reynolds  shear  stress. 
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The  relative  importance  of  the  third  term,  and  to  a  much  lesser  extent  the  fourth,  is  based 
somewhat  on  injector  parameters.  Steeper  injection  angles  or  higher  injection  velocities  will 
increase  both  terms. 

The  multiple  overheat  data  validates  the  T;  =  0  assumption  made  in  the  present 
work.  It  also  indicates  the  single  overheat  data  to  be  accurate  within  approximately  3%  of 
the  multiple  overheat  data.  For  flowfields  with  small  total  temperature  fluctuations,  single 
overheat  data  provides  excellent  accuracy  (within  3%). 

6.2  Recommendations  for  Further  Study 

The  importance  of  compressible  turbulence  in  supersonic  injection  cannot  be 
overemphasized.  More  detailed  multiple  overheat  analysis  is  needed.  In  fact,  it  may  be 
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possible  to  devise  a  correlation  between  single  and  multiple  overheat  hot-wire  data  to  extend 

the  usefulness  of  the  single  overheat  method. 

Innumerable,  useful  investigations  can  be  undertaken  with  the  AFIT  Mach  3  wind 
tunnel  and  low  angle  injection  model.  Heated  injectant  can  be  used  to  determine  the  effect 
of  differences  in  total  temperature  on  turbulence  and  mixing  (this  work  is  underway  as  of 
this  writing).  Low  molecular  weight  injectants  can  be  used  to  simulate  injection  in  a 
hydrogen-fueled  SCRAMjet  engine.  High  molecular  weight  injectant  can  simulate  the 
conditions  found  in  a  chemical  oxygen-iodine  laser  (currently  under  development  at  Phillips 
Lab  in  Albuquerque,  NM). 

Better  investigative  techniques  can  be  used.  Seeded  flows  with  laser  doppler 
velocimetry  (LDV)  may  yield  new  insights.  Use  of  triple-wire  probes  can  reduce  tunnel  time 
and  increase  fidelity  of  the  flow  components.  In  particular,  the  triple-wire  would  enable  the 
measurement  of  the  one  element  of  the  shear  stress  tensor  not  measured  here  ( 


Tunnel  run  time  can  be  increased  by  addition  of  extra  vacuum  tanks.  While  doubling 
the  vacuum  storage  would  double  the  mean  time  between  tunnel  runs,  the  useful  tunnel  time 
would  more  than  double  because  the  fixed  costs  of  starting  the  tunnel,  allowing  it  to  settle, 
and  triggering  the  DAS  would  only  have  to  be  paid  once. 

Finally,  to  enhance  further  study  and  applicability  of  future  data,  the  shock  structure 
must  be  more  fully  understood.  This  can  be  accomplished  by  LDV,  vertical  schlieren  or 
shadowgraph  photography,  or  by  traversing  the  mean  flow  probes  in  horizontal,  instead  of 

vertical,  sections. 
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